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ABSTRACT 

The recombination of hydrogen and helium at z ~ 1000-7000 gives unavoidable distor- 
tions to the Cosmic Microwave Background (CMB) spectrum. We present a detailed 
calculation of the line intensities arising from the Ly a (2p-ls) and two-photon (2s-ls) 
transitions for the recombination of hydrogen, as well as the corresponding lines from 
helium. We give an approximate formula for the strength of the main recombination 
line distortion on the CMB in different cosmologies, this peak occurring at about 
170 /im. We also find a previously undescribed long wavelength peak (which we call 
the pre-recombination peak) from the lines of the 2p-ls transitions, which are formed 
before significant recombination of the corresponding atoms occurred. Detailed calcu- 
lations of the two-photon emission line shapes are presented here for the first time. 
The frequencies of the photons emitted from the two-photon transition have a wide 
spectrum and this causes the location of the peak of the two-photon line of hydrogen 
to be located almost at the same wavelength as the main Ly a peak. The helium lines 
also give distortions at similar wavelengths, so that the combined distortion has a com- 
plex shape. The detection of this distortion would provide direct supporting evidence 
that the Universe was indeed once a plasma. Moreover, the distortions are a sensitive 
probe of physics during the time of recombination. Although the spectral distortion 
is overwhelmed by dust emission from the Galaxy, and is maximum at wavelengths 
roughly where the cosmic far-infrared background peaks, it may be able to tailor an 
experiment to detect its non-trivial shape. 

Key words: lines: formation - cosmology: cosmic microwave background - cosmol- 
ogy: early universe - cosmology: theory - atomic processes - infrared: general. 



1 INTRODUCTION 

Physical processes in the plasma of the hot early Universe 
thermalize the radiation content, and this redshifts to be- 
co me the observed Cosm ic Microwave Background (CMB; 
see IScott fe Smootl 120041 and references therein) . Besides 
the photons from the radiation background, there were some 
extra photons produced from the transitions when the elec- 
trons cascaded down to the ground state after they recom- 
bined with the ionized atoms. The transition from a plasma 
to mainly neutral gas occurred because as the Universe ex- 
panded the background temperature dropped, allowing the 
ions to hold onto their electrons. The photons created in this 
process give a distortion to the nearly perfect blackbody 
CMB spectrum. Since recombination happens at redshift 
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z ~ 1000, then Lya is observed at ~ 100/im today. There 
is approximately one of these photons per baryon, which 
should be compared with the ~ 10 9 photons per baryon in 
the entire CMB. However, the recombination photons are 
superimposed on the Wien part of the CMB spectrum, and 
so make a potentially measurable distortion. 

From the Far-Infrared Absolute Spe c tropho tome- 
ter (FIRAS) mea surements, iFixsen et alJ il996T) and 
iMather et all J199ST) showed that the CMB is well modelled 
by a 2.725±0.001 K blackbody, and that any deviations from 
this spectrum around the peak are less than 50 parts per mil- 
lion of the peak brightness. Constraints on smooth functions, 
such as fj,- or ^-distortions are similarly very stringent. How- 
ever, there are much weaker constraints on narrower features 
in the CMB spectr um. Moreover with in the last decade it 
has been discovered iPueet et al.l lll996[) that there is a Cos- 
mic Infrared Background (CIB; see Mauser fc Dwekl l200ll . 
and references therein), which peaks at 100-200//m and is 
mainly composed of luminous infrared galaxies at moderate 
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redshifts. The existence of this background makes it more 
challenging to measure the recombination distortions than 
would have been the case if one imagined them only as be- 
ing distortions to Wien tail of the CMB. However, as we 
shall see, the shape of the recombination line distortion is 
expected to be much narrower than that of the CIB, and 
hence the signal may be detectable in a future experiment 
designed to measure the CIB spectrum in detail. 

The first published calculations of the line distortions 
occur in the semi nal pape r s on t he cosmological recombi- 
nation process by IPeeblesI 1^)68) and Zel'dovich, Kurt & 
Sunyaev (1968). One of the main motivations for study- 
ing the recombination process was to answer the question: 
'Where are the Ly a line photons from the recombination in 
the Universe?' (as reported in Rubino-Martin, Hernandez- 
Monteagudo & Sunyaev 2005). In fact these studies found 
that for hydrogen recombination (in a cosmology which is 
somewhat different than the model favoured today) there are 
more photons created through the two- photon 2s-ls tr ansi- 
tion than from t h e Ly a transition. Both lPeeblesI dl968f) and 
IZel'dovich et alJ (Il968l) plot the distortion of the CMB tail 
caused by these line photons, but give no detail about the 
line shapes. Other authors have included some calculation 
or discussion of the line disto rtions as part of other recom- 
bination related st udies, e.g. iBosc han fc Biltzingerj lll998l) . 
and most recently ISwitzer & Hh'ateril2005l) . However, the 
explicit line shapes have never before been presented, and 
the helium lines have also been neglected so far. The only 
numerical study to show th e hydrogen lines in any detail i s 
a short conference report bv lDeH'Antonio fc Rvbickil p993), 
meant as a preliminary version of a more full study which 
never appeared. Although their calculation appears to have 
been substantially correct, unfortunately in the one plot 
they show of the distortions (their fig. 2) it is difficult to 
tell precisely which effects are real and which might be nu- 
merical. 

Some of the recombination line distortions from higher 
energ y levels, n > 2, have als o been c alculated iDubrovichl 
19751: iLvubarskv fc Sunyaevl Il983l: iFahr fc Lochl Il99lf 
Burdvuzha fc Chekmezovl Il994t iDeir Antonio fc Rybickl 
19931: iDubrovich fc Stolvarovl Il995l Il997l: iBurginl l2003t 
Kholupenko. Ivanchik fc Varshalovichl I2005F) . However, 
these high n lines are extremely weak compared with the 
CMB (below the 10" 6 level), while the Lya line is well 
above the CMB in the Wien region of the spectrum. 

As trumpeted by many authors, we are now entering 
into the era of precision cosmology. Hence one might imag- 
ine that future delicate experiments may be able to mea- 
sure these line distortions. Since the lines are formed by the 
photons emitted in each transitions of the electrons, they 
are strongly dependent on the rate of recombination of the 
atoms. The distortion lines may thus be a more sensitive 
probe of recombination era physics than the ionization frac- 
tion x e , and the related visibility function which affects the 
CMB anisotropies. This is because a lot of energy must be 
injected in ord er for any physical process t o change x e sub- 
stantially fe.g. iPeebles. Seager fc H u 2000). In general that 
energy will go into spectral distortions, including boosting 
the recombination lines. 

This also means that a detailed understanding of the 
physics of recombination is crucial for calculating the dis- 
tortion. The basic physical picture for cosmological recom- 



bination ha s not changed si n ce the early work of iPeeblea 
lll968ft and IZel'dovich et alJ dl968f) . However, there have 
been several refinements introduced since then, motivated by 
the increased emphasis on obtaining an accurate recombina- 
tion history as part of the calculation of CMB anisotropies. 
Seager, Sasselov & Scott (1999,2000) presented a detailed 
calculation of the whole recombination process, with no 
assumption of equilibrium among the energy levels. This 
multi-level computation involves 300 levels for both hy- 
drogen and helium, and gives us the currently most accu- 
rate p icture of the recom bination history. In the context 
of the ISeager et alJ (I200CT) recombination calculation, and 
with the well-developed set of cosmological parameters pro- 
vided by Wilkinson Microwave Anisotropy Probe (WMAP; 
ISpergel et al l2003h and other CMB experiments, it seems an 
appropriate time to calculate the distortion lines to higher 
accuracy in order to investigate whether they could be de- 
tected and whether their detection might be cosmologically 
useful. 

The aim of this paper is to calculate the line distortions 
on the CMB from the 2p-ls and 2s-ls transitions of H and 
the corresponding lines of He (i.e. the 2 p-1 s and 2 s-l 1 s 
transitions of Hel, and the 2p-ls and 2s-ls transitions of 
Hell) during recombination, using the standard cosmologi- 
cal parameters and recombination history. In Section [5] we 
will describe the model we used in the numerical calcula- 
tion and give the equations used to calculate the spectral 
lines. In Section [3] we will present our results and discuss 
the detailed physics of the locations and shapes of the spec- 
tral lines. An approximate formula for the magnitude of the 
distortion in different cosmologies will also be given. Other 
possible modifications of the spectral lines and their poten- 
tial detectability will be discussed in Section 2] And finally, 
we will present our conclusions in the last section. 

2 BASIC THEORY 
2.1 Model 

Instead of adopting a full multi-level code, we use a sim- 
ple 3-level model atom here. For single-electron atoms (i.e. 
HI and Hell), we consider only the ground state, the first 
excited state and the continuum. For the 2-electron atom 
(Hel), we consider the corresponding levels among singlet 
states. In general, the upper level states are considered to be 
in thermal equilibrium with the first excited state. Case B re- 
combination is adopted here, which means that we ignore re- 
combinations and photo-ionizations directly to ground state. 
This is because the photons emitted from direct recombi- 
nations to the ground state wil l almost immediately reion- 
ize a nearby neutral H atom JPeeblesl fl96^ : ISeaeer et all 
2000). We also include the two-photon rate from 2s to 
the ground st ate for all three atoms, with rates: A^_ 1k = 
8.229 063 s" 1 JGoldmanlll98l ISantos. Parente fc Indelicate! 
Il998ft : A« 1 e s I _ lla = 51.02 s- ^Derevianko fc Johnson! 11997^ 
although it makes no noticeable difference to the calcula- 
tion if one uses the older value of 51.3 s _1 from Drake, 
Victor fc Dalgamrno (1969) ; and A? a e l\ a = 526.532s -1 
jLipeles. Novick fc TolkllT965l lGoldmarJll989ll . This 3-level 
atom model is similar to the one used in the program recfast, 
with the main difference being that here we do not assume 
that the rate of change of the first excited state 712 is zero. 
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The rate equations for the 3 atoms are similar, and so 
we will just state the hydrogen case as an example: 
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Here the values of n; are the number density of the ith state, 
where n e and n p are the number density of electrons and 
protons respectively. AR^-j is the net bound-bound rate 
between state i and j and the detailed form of AR^ p _ lB and 
Ai?JJ_ ls will be discussed in the next subsection. H(z) = 
a/ a is the Hubble factor, 
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where Q represents the fraction of the critical density in 
matter, curvature or cosmological constant, and the Hubble 
parameter today Ho = lOO/i km s~ Mpc ~ . Finally qh i s 
the Case B recombination coefficient from iHummeJ lll994t) . 
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1 + ct d 

which is fitted by IPeauignot. Petitiean fc I (|l99lfl , with a = 
4.309, b = -0.6166, c = 0.6703, d = 0.5300 and t = 
Tm/10 4 K, while /3h is the photo-ionization coefficient: 
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where fee is Boltzmann's constant, h p is Planck's constant, 
m c is the mass of electron, Tm is the matter temperature 
and i^2s-c is the frequency of the energy difference between 
state 2s and the continuum. For the rate of change of Tm, 
we include only th e Compton and adiabatic cooling terms 
JSeaeer et alJl2000D . i.e. 



8a T U 



q . x dT M _ 

dz 3H(z)m e c n e + nu + nnc 
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where Tr is the radiation temperature, c is the speed of 
light, U — chrTr, or is the radiation constant and crx is the 
Thompson scattering cross-section. 

We use the Bader-Deuflhard semi-impli cit numerical in- 
tegration scheme (see Section 16.6 in lPress et aljfl99^ to 
solve the above rate equations. All the numerical results are 
made using the ACDM model with parameters: Qb = 0.046; 
n M = 0.3; fl A = 0.7; Ok = 0; Y v = 0.2 4; T = 2.725 K 
and h = 0.7 (see e.g. ISpergel et alll2003h . Here Y p is the 
primordial He abundance and To the present background 
temperature. 



2.2 Spectral distortions 

We want to calculate the specific line intensity I VQ (z — 0) 
(i.e. energy per unit time per unit area per unit frequency 
per unit solid angle, measured in Wm _2 Hz _1 sr _1 ) observed 
at the present epoch, z — 0. The detailed calculation of 
Iv (z — 0) for the Lya transition and the two-photon tran- 
sition in hydrogen a re presented as exam ples (the notation 
follows Section 2.5 in IPadmanabhanl 19931) . A similar deriva- 
tion holds for the corresponding transitions in helium. To 
perform these calculations we first consider the emissivity 
jv{z) (energy per unit time per unit volume per unit fre- 
quency, measured in Wm _3 Hz _1 ) of photons due to the 
transition of electrons between the 2p and Is states at red- 
shift z: 



j u {z) = h p vAR% p _ le ,{z)4>[v(z)\, 



(9) 



where 4>(y) is the frequency distribution of the emitted pho- 
tons from the emission process and AJ?2p-is is the net rate 
of photon production between the 2p and Is levels, i.e. 



AR% p _ ls = pi2 (njpifei - ni R12) 



(10) 



Here is the number density of hydrogen atoms having 
electrons in state i, the upward and downward transition 
rates are 



and 



R12 = B12J, 

R21 = (.421+B21J) 



(11) 
(12) 



with A21, B\2 and B21 being the Einstein coef ficients an d 
p\2 the Sobolev escape probability fsee lSeaeer et alJ l2000). 
which accounts for the redshifting of the Lya photons due 
to the expansion of the Universe. As nf 2> ri§ pl can be 
expressed in the following form: 



P12 



, with 
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We approximate the background radiation field J as a per- 
fect blackbo dy spectrum by ig noring the line profile of the 
emission fsee lSeaeer et al]l2000T) . We also neglect secondary 
distortions to the radiation field (but see the discussion in 
Section I4.1H . These secondary distortions come from pho- 
tons emitted earlier in time, during recombination of H or 
He, primarily the line transitions described in this paper. 
Assuming a blackbody we have 

rfrp \ 2hpV a 

JyJ-Mj = 5 — 



exp 



/ h p u a 



(15) 

where u a = c/121.5682 nm= 2.466 x 10 15 Hz and corre- 
sponds to the energy difference between states 2p and Is, 
while the frequency of the emitted photons is equal to v a . 
Therefore, we can set <f)[v(z)] = 5[v(z) — v a ], i.e. a delta 
function centred on v a , so that 



jl ya {z) = h p vAR^ p _ 



(z)8[u( 



(16) 



The increment to the intensity coming from time interval dt 
at redshift z is 



dl„{z) 



In 



jvdt, 



(17) 



which redshifts to give 
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dl„ (z = 0) 



4tt (1 + z) 



-dt. 



We assume that the emitted photons propagate freely until 
the present time. Integration over frequency then gives 



with 
l + z a 



0) 
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Chm 



-dt 



(1 + z) 3 
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(18) 
(19) 



f ' 



using 

i/(z) = fo(l + z) 



and 



dt 

dz 



H(z)(l + z)- 

Equation 11911 is the basic equation for determining the Ly a 
line distortion, using AR^ p _ ls (z) from the 3-level atom cal- 
culation. 

For the two-photon emission between the 2s and Is levels, 
the emissivity at each redshift is 



j u (z) = hpvAR 2s -i 3 (z)cf>lv(z)], 



(20) 



and the calculation is slightly more complicated, since for 
4>{v) we need the frequency spectrum of the emission pho- 
tons of the 2s-ls trans i tion of H dSpitzer fc Greensteinl 
Il95ll: iMartinis fc Stoiid l2000l) as shown in Fig. □ Here 
Ai?^ s _ ls is the net rate of photon production for the 2s- 
ls transition, i.e. 



Ai? 2s _ ls = A H (n 2s - rite p a/ B M j . 
Therefore, using equation 11811 . we have 



iVA* = o) 



ch v uo 
4ir 



ARg. 



*)#*(! + *)] 



H{z)(l + z) 



(21) 



dz. (22) 



We use the sim ple trapezoidal rule (see Section 4.1 in 
iPress et al.l Il992h to integrate equation 1221 numerically 
from 2 = to the time when AR is sufficiently small that 
the integrand can be neglected. 



3 RESULTS 

Each of the line distortions are shown separately in Fig. [5] 
and summed for each species in Fig. [3] The shape of the 
lines from H, Hel and Hell are fairly similar. There are two 
distinct peaks to the 2p-ls emission lines. We refer to the 
one located at longer wavelength as the 'pre-recombination 
peak', since the corresponding atoms had hardly started to 
recombine during that time. The physics of the formation 
of this peak will be discussed in detail in section Ti. 1.1 1 The 
second (shorter wavelength) peak is the main recombination 
peak, which was formed when the atoms recombined. While 
the longer wavelength peak actually contains almost an or- 
der of magnitude more flux, it makes a much lower relative 
distortion to the CMB. The ratio of the total distortion to 
the CMB intensity is shown in Fig. 0] It is approximately 
one for the main recombination peak, but ~ 10~ 4 for the 
pre-recombination peak. 

In Fig. |21 we plot the lines from H and Hel together with 




Figure 1. The normalized emiss ion spectrum for the two- 
photon process (2s— Is) of hydrogen iSpitzer fc Greensteinlll95lt 
IMartinis fc Stonal200(]() . The top panel shows </>[j/(A)] vs A, while 
the bottom panel shows <j>{y) vs y, where v = yu a . Note that the 
spectrum is symmetric in v about v a /2, but the A spectrum is 
very asymmetric, being zero below A a , and having a tail extend- 
ing to high A. 



H line 
Hel line 
Hell line 




Figure 2. The line intensity \oI\ from the net Lya emission 
of H (thick solid), the two-photon emission (2s-ls) of H with the 
spectrum <j>{v) (thin solid), the 2 1 p— l x s emission of Hel (thick 
dashed), the 2 1 s-l 1 s two-photon emission of Hel (thin dashed), 
the 2p-ls emission of Hell (thick dotted) and the 2s— Is two- 
photon emission of Hell (thin dotted). 



the CMB and an estimate of the CIB. We can see that the 
lines which make the most significant distortion to the CMB 
are the Lya line and the 2 1 p-l 1 s line of Hel, and that these 
lines form a non-trivial shape for the overall distortion. The 
sum of all the spectral lines and the CMB is shown in Fig. |S| 
Note that these lines will also exist in the presence of the 
CIB - but the sha pe of this background is currently qui te 
poorly determined jFixsen et alJll998l : lHauser et aTlll998l) . 

We now discuss details of the physics behind the shapes 
of each of the main recombination lines. 



Spectral distortions to the CMB 5 




Figure 3. The line intensity \qI\ from the sum of the net Lya 
emission and two-photon emission (Is— 2s) of H (thick solid), the 
sum of the 2 1 p— l^s emission and 2 1 s-l 1 s two-photon emission 
of Hel (thick dashed), and the sum of the 2p-ls emission and 
2s-ls two-photon emission of Hell (thick dotted), together with 
the backgrou nd spectra: CMB ( long-dashed); and estimated CIB 
(dot-dashed: iFixsen et aT1ll99Sft . The sum of all the above emis- 
sion lines of H and He plus the CMB is also shown (thin solid). 




Figure 5. Comparison of the net 2p— Is (solid) and 2s— Is (dashed) 
transition rates of H. The Lya redshifting process dominates 
during the start of recombination, while the 2-photon process is 
higher during most of the time that recombination is occurring. It 
turns out that in the standard ACDM model about equal numbers 
of hydrogen atoms recombine through each process, with slightly 
over half the hydrogen in the Universe recombining through the 
2-photon process. 




Figure 4. The ratio of the total line distortion to the CMB inten- 
sity is plotted. The ratio is larger than 1 (i.e. the intensity of the 
distortion line is larger than that of the CMB) when Ao ~ 170 fim 
which is just where the main Lya line peaks. 

3.1 Lines from the recombination of hydrogen 

During recombination, the Lyman lines are optically thick, 
which means that photons emitted from the transition to 
n = l are instantly reabsorbed. However, some of the emit- 
ted photons redshift out of the line due to the expansion of 
the Universe and this makes the Lya transition one of the 
possible ways for electrons to cascade down to the ground 
state. The other path for electrons going from n = 2 to 
n = 1 is the two-photon transition between 2s and Is. Fig. 
shows the net photon emission rate of the Lya and two- 
photon transitions as a function of redshift for the stan- 
dard ACDM model. The two-photon rate dominates at low 
redshift, where the bulk of the recombinations occur. This 
means that there are more photons emitted through the 
two-photon emission process (54% of the total number of 



photons created during recombination of H) than through 
the Lya redshif t ing p rocess. This conclusion agrees with 
IZel'dovich et al.l l|l96ftl) - although of course the balance 
depen ds on the cosmological parameters (see ISeaeer et alJ 
2000) and for today's best fit cosmology the two processes 
are almost equal. Despite this fact, the overall strength of 
the two-photon emission lines are weaker because the pho- 
tons are not produced with a single frequency, but with a 
wide spectrum ranging from to v a . The location of the two- 
photon peak (see Fig.[5J is also somewhat unexpected, since 
it is almost at the same wavelength as the Ly a recombina- 
tion peak, rather than at twice the wavelength. The reason 
for this will be discussed in section l3.1.2l 

We should also note that the tiny dip in our curves for 
the long-wavelength tail of the pre-recombination peak (see 
Fig-HJ is due to a numerical error, when the number density 
of the ground state is very small. This can also be seen in 
the pre-recombination peak for Hell. 



3.1.1 The pre-recombination emission peak 

The highest Ly a peak (shown in Fig. |2J is formed before 
the recombination of H has already started, approximately 
at z > 2000. During that time the emission of Ly a photons 
is controlled by the bound-bound Ly a rate from n = 2 (i.e. 
the TI2-R21 term in equation 1101 ) and the photo- ionization 
rate (712 «h). From Fig. [(J we can see that at early times the 
bound-bound Lya rate is larger than the photo- ionization 
rate. This indicates that when an electron recombines to 
the n = 2 state, it is more likely to go down to the ground 
state by emission of a Ly a photon than to get ionized. The 
excess Ly alpha photons are not reabsorbed by ground state 
H, but are redshifted out of the absorption frequency due 
to the expansion of the Universe; they escape freely and 
form the pre-recombination emission line. Note that there is 
very little net recombination of H, since the huge reservoir 
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Figure 6. The top panel shows the bound-bound Ly a rate n^Rix 
and the photo-ionizing rate n2«H for n=2. The lower panel shows 
the fraction of ground state H atoms ni/njj, and also the ioniza- 
tion fraction x e . 



of > 13.6 eV CMB photons keeps photo-ionizing the ground 
state H atoms (see Fig. I12> . 

We now turn to a more detailed explanation of the pre- 
recombination emission peak. The bound-bound Lya rate 
from n = 2 is initially approximately constant, as it is domi- 
nated by the spontaneous de-excitation rate (the A21 term in 
equation l|21flL At the same time the photo- ionization rate 
is always decreasing as redshift decreases, since the num- 
ber of high energy photons keeps decreasing with the ex- 
pansion of the Universe. Therefore, with a constant bound- 
bound Lya rate and the decreasing photo-ionization rate, 
the emission of Lya photons rises. The peak of this pre- 
recombination line of H occurs at around z = 3000, by which 
time only a very tiny amount of ground state H atoms have 
formed (n\/nn < 10 -7 , see Fig. These ground state H 
atoms build up until they can reabsorb the Lya photons 
and this lowers the bound-bound Lya rate. The decrease 
of the bound-bound Ly a rate is represented in the Sobolev 
escape probability pi2 in equation 1131 . At high redshift, 
P12 is 1 and there is no trapping of Lya photons. When H 
starts to recombine, the optical depth r s increases and the 
Ly a photons can be reabsorbed by even very small amounts 
of neutral H. For r s 3> 1, we can approximate P12 — l/r s 
and P12 oc H{z)/ni. Because of the increase in the number 
density of the ground state and the decrease in H(z), the 
pre-recombination line decreases. One can therefore think 
of the 'pre-recombination peak' as arising from direct Lya 
transitions, before enough neutral H has built up to make 
the Universe optically thick for Lyman photons. This process 
occurs because the spontaneous emission rate (A21 term) is 
faster than the photoionization rate for n — 2; it increases 
as the Universe expands, due to the weakening CMB black- 
body radiation, and is quenched as the fraction of atoms in 
the n = 1 level grows. The shorter wavelength peak, on the 
other hand, comes from the process of redshifting out of the 
Ly a line during the bulk of the recom bination epoch. 

By using the RECFAST program JSeager et alJll99Sft . 
we can generate the main Ly a recombination peak and also 
the two-photon emission spectrum, by simply adding a few 
lines into the code. However, the pre-recombination peak 



cannot be generated from RECFAST, since there the rate of 
change of the number density of the first excited state 712 
is assumed to be negligible and is related to n\ via thermal 
equilibrium. Moreover, in the effective 3-level formalism, the 
Lya line is assumed to be optically thick throughout the 
whole recombination process of H (in order to reduce the 
calculation into a single ODE), which is not valid at the 
beginning of the recombination process. Hence, one needs 
to follow the rate equations of both states (i.e. n = 1 and 
n = 2) to generate the full Lya emission spectrum. The 
pre-recombination peak of H was mentioned and p lotted in 
the earlier work of iDell' Antonio fc Rvb icki ( 1993j) as well, 
although they did not describe it in any detail. 

Another way to understand the line formation mech- 
anism is to ask how many photons are made in each pro- 
cess per atom. We find that for the main Lya peak there 
are approximately 0.47 photons per hydrogen atom (in the 
standard cosmology). During the recombination epoch, net 
photons for the n — 2 to n = 1 transitions are only made 
when atoms terminate at the ground state. Hence we expect 
exactly one n — 2 to n — 1 photon for each atom, split be- 
tween the Lya redshifting and 2-photon processes (and the 
latter splits the energy into two photons, so there are 1.06 of 
these photons per atom). For the 'pre-recombination peak', 
on the other hand, the atoms give a Ly a photon when they 
reach n = 1, but they then absorb a CMB continuum pho- 
ton to get back to higher n or become ionized. The number 
of times an atom cycles through this process depends on 
the ratio of the relevant rates. If we take the rate per unit 
volume from Fig. |S] and divide by the number density of hy- 
drogen atoms at z ~ 3000 then we get a rate which is about 
an order of magnitude larger than the Hubble parameter 
at that time. Hence we expect about 10 'pre-recombination 
peak' photons per hydrogen atom. A numerical calculation 
gives the more precise value of 8.11. 

3.1.2 The two-photon emission lines 

Surprisingly, the location of the peak of the line intensity of 
the 2s-ls transition is almost the same as that of the Lya 
transition, as shown in Fig.|2] while one might have expected 
it to differ by a factor of 2. In order to understand this effect, 
we rewrite the equation (1221 in the following way: 

/>oo 

O* = °) = / = 0; *0 dz\ (23) 

Jo 

where 4>(z') — uo<p(u'), and 

^(, = 0^-^ = 0;^)) = ^^^, (24) 
with 

l + z>^-. 

Equation 1'2 H gives the redshifted flux (measured now at 
z — 0) of a single frequency v' coming from redshift z' and 
corresponding to the redshifted frequency fg. 

We first calculate the line intensity of the two-photon 
emission with a simple approximation: a delta function spec- 
trum S(v — v a /2), where v a /2 is the frequency corresponding 
to the peak of the two-photon emission spectrum <f>(v). Fig.0 
shows the intensity spectrum of two-photon emission using 
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v„(Hz) 



Figure 7. The line intensity of the 2s-ls transition (two-photon 
emission) Iv (z = 0) as a function of redshifted frequency uo 
for three different assumptions: the correct frequency spectrum 
of two-photon emission (solid); the delta function approximation 
8(u — u a /2) (dashed); and the fiat spectrum approximation (dot- 
ted). 




l+z' 



Figure 8. The top panel shows the redshifted flux from single 
emission frequency li {z = 0;z') plotted against the redshift of 
emission, l + z' . The bottom panel shows the frequency spectrum 
of two-photon emission 4>[i/(z')] plotted against z' for three red- 
shifted frequencies: u = 10 12 Hz; 1.6 X 10 12 Hz; and 5 X 10 12 Hz. 

a delta frequency spectrum 8(u — v a /2) compared with the 
two-photon emission using the correct spectrum <f>(v) . We 
can see that there is a significant shift in the line centre com- 
pared with the 5-function case. Where does this shift come 
from? 

We know that the frequencies of emitted photons are 
within the range of to u a at the time of emission. For 
a fixed redshifted frequency v§ now, we can calculate the 
range of emission redshifts contributing to uo (referred to 
as the 'contribution period' from now on), which is repre- 
sented by (j)'(z') or 4>(y 1 ') . In Fig. [SJ we show the spectral 
distribution cj>[v'{z')] as a function of redshift z for specific 
values of uo- For example, if we take vq — 5 x 10 12 Hz, 
then photons emitted between 1 + z = 1 (i.e. v = v§) and 
~ 500 (y — u a ) will give contributions to vq. The smaller the 
redshifted frequency vq, the wider the contribution period. 



We might expect that the line intensity of this two-photon 
emission will be larger if the contribution period is longer, as 
there are more redshifted photons propagating from earlier 
times. However, this is not the case, because the rate of two- 
photon emission R,2-y also varies with time, and is sharply 
peaked at z ~ 1300-1400. Hence (z — 0; z') is also sharply 
peaked at z ~ 1300-1400. In Fig. |H| the redshifted flux in- 
tegrand lt (z — 0, z) and the emission spectrum are 
plotted on the same redshift scale. For vq = 5 x 10 12 Hz (low- 
est panel), we can see that the contribution period covers a 
redshift range when J* (z = 0, z) and i?2 7 are small in value. 
The contribution period widens with decreasing i/q and cov- 
ers more of the redshift range when two-photon emission is 
high. Therefore, the flux I vo {z = 0) is expected to increase 
with decreasing vq until the contribution period extends to 
the redshifts at which the two-photon emission peaks. As 
i/q gets even smaller (e.g. vq — 10 12 Hz), then the contribu- 
tion period becomes larger than the redshift range for two- 
photon emission and hence only lower energy photons can 
be redshifted to that redshifted frequency. As a result, the 
flux I VQ (z — 0) starts to decrease, and so we have a peak. 
The flux peaks at vq ~ 10 12 Hz when we use the 5-function 
approximation. However, from Fig. |H| we can see that the 
contribution period for uo — 10 12 Hz is much greater than 
that of the two-photon emission period, and therefore this is 
not the location of peak. Based on the argument presented 
above, we expect the peak to be at around 1.6 x 10 12 Hz, or 
200 Atm. 

The basic mathematical point is that 4>{y) is extremely 
poorly represented by a 5-function. Since the spectrum <^>(^) 
is quite broad, it can be better approximated as a uniform 
distribution than as a 5-function. Another crude approxima- 
tion would be to assume a flat spectrum for 4>{v) in Fig. 
Fig.[7]compares the intensity I VQ (z — 0) found using the cor- 
rect form for 4>{v) with the 5-function and flat spectrum ap- 
proximations. This shows that the flat spectrum gives quali- 
tatively the same results as the correct form of the spectrum, 
and that the peak occurs fairly close to that of Lya, but is 
much broader. The same general arguments apply to the 
two-photon lines of Hel and Hell (as we discuss in Section 
3.2). 

3.1.3 Dependence of Qm and Qb 

The largest distortion on the CMB is from the shorter wave- 
length recombination peak of the hydrogen Lya line (see 
Fig. It may therefore be useful to estimate the peak of 
this line's intensity as a function of the cosmological pa- 
rameters. The relevant parameters are the matter density 
(oc Hm/i 2 ) and the baryon density (oc Qsh 2 ). This is be- 
cause Q,mh 2 affects the expansion rate, while Qsh 2 is related 
to the number density of hydrogen. No other combinations 
of cosmological parameters have a significant impact on the 
physics of recombination. 

We can crudely understand the scalings of these pa- 
rameters through the following argument. Regardless of 
the escape probability pi2, the remaining part of the rate 
( n 2p^2i — rii -R12) is roughly proportional to n^ 1 oc Q,Bh 2 (l — 
x e ). The escape probability P12 can be approximated as 1 at 
the beginning of recombination (r s <C 1) and l/r s during the 
bulk of the recombination process (with r s 3> 1). Note that 
r s oc H(z)/nf oc {n M h 2 ) 1/2 [D. B h 2 {l - Therefore, 
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A_R 2 p-ls OC 

and thus 

AR 



{n M h 2 ) [n B h 2 (i~x e )] 
(n M h 2 ) 1/2 [n B h 2 (i-x e 



for r s < 1 
1° for r s > 1, 



(25) 



H{z) 



(n M h 2 )- 1/2 {n B h 2 (i- 

(Q M h 2 )°[Q B h 2 (l - x e 



for r s < 1 
for r s > 1. 



(26) 



From this rough scaling argument, we may expect that the 
SIm dependence of the peak of the Lya line is an approxi- 
mate power law with index between —1/2 and 0, while for 
Q B the corresponding power-law index is expected to lie be- 
tween and 1. The dependence of Hm is actually more com- 
plicated when one allows for a wider range of values (see 
iDell' Antonio fc Rvbickil ll993), The above estimation just 
gives a rough physical idea of the power of the dependence. 

A more complete numerical estimate of the peak of the 
recombination Lya distortion is: 



(Aq/a ) 



peak 



: 8.5 x 10" 



Q B h 2 V 
0.0224 J 



0.147 



Wm sr 



(27) 




Figure 9. Comparison of the net 2 1 p— l^^s (solid) and 2 1 s—l 1 s two- 
photon (dashed) transition rates of Hel. The two-photon rate is 
sub-dominant through most of the Hel recombination epoch, and 
hence, unlike for hydrogen, most helium atoms did not recombine 
through the two-photon process. 



where we have normalized to the parameters of the currently 
favoured cosmological model. The peak occurs at 

A ~170Atm (28) 

for all reasonable variants of the standard cosmology. 

3.2 Lines from the recombination of helium (Hel 
and Hell) 

We compute the recombination of Hell and Hel in the same 
way as for hydrogen. For the two-electron atom Hel, we ig- 
nore all the forbidden transitions between singlet and triplet 
states due to the low p opulation of the triplet states (see 
ISeaeer et alJll999ll2000l) . The 2 1 p-l 1 s transitions of Hel are 
optically thick, the same situation as for H. This makes the 
electrons take longer to reach the ground state and causes 
the recombination of Hel to be slower than Saha equilibrium. 
However, unlike for H, and despite the optically thick 2 1 p- 
l 1 s transition line, the 2 1 p-l 1 s rate dominates, as shown in 
Fig. |U] For Hell, due to the fast two-photon transition rate 
(see Fig. I1UI . there is no 'bottleneck' at the n — 2 level in 
the recombination process. Hence Hell recombination can 
be we ll approximated b y using the Saha equilibrium for- 
mula JSeaeer et alJl200fJ) . 

We can see the effect of the above differences in recom- 
bination history on the lines: the width of the recombination 
peak of both H and Hel is larger than that of Hell. Overall, 
the spectral lines of Hell are of much lower amplitude than 
those of H (see Fig.|2J with the distortion to the CMB about 
an order of magnitude smaller. 

The peaks of the line distortions from H and Hell are 
located at nearly the same wavelengths. For hydrogenic ions 
the ls-2p energy (and all the others) scales as Z 2 , where Z 
is the atomic number. Hence for Hell recombination takes 
place at z~ 6000 rather than the z ~ 1500 for hydrogen. 
Hence the line distortion from the 2p-ls transition of Hell 
redshifts down to about 200 fim, just like Ly a. 

The two-photon frequency spectrum of Hell is the 
same as for H, since they both are single-electron 



atoms llTung et al.lll984l) . However, it is complicated to cal- 
culate the two-photon frequency spectrum of Hel very ac- 
curately, sin ce there is no exact wave-function for the state 
of the atom. iDrake et alJ lll969T) used a variational method 
to calculate the two-photon frequency sp ectrum of Hel with 
values given up to 3 significant figures. Drake ( 1986) pre- 
sented another calculation, giving one more digit of pre- 
cision, and making the spectrum smoother, as shown in 
Fig. 1111 These two calculations differ by only about 1%, 
which makes negligible change to the two-photon Hel spec- 
tral line. 

All of the H and He lines (for n = 2 to n = 1) are 
presented in Fig. [5] and the sum is shown as a fractional 
distortion to the CMB spectrum in Fig. [I] We find that 
in the standard cosmological model, for Hel recombination, 
there are about 0.67 photons created per helium atom in 
the 'main' 2 1 p-l 1 s peak, 0.70 per helium atom in the 'pre- 
recombination peak', and 0.66 in the two-photon process. 
The numbers for Hell recombination are 0.62, 0.76 and 6.85 
for these three processes, respectively. 



4 DISCUSSION 

4.1 Modifications in the recombination calculation 

There are several possible improvements that we could make 
to the line distortion calculation. However, as we will discuss 
below, we do not believe that any of them will make a sub- 
stantial difference to the amplitudes of the lines. 

In order to calculate the distortion lines to higher accu- 
racy, we should use the multi-level model without any ther- 
mal equilibrium assumption among the bound states. And 
we also need to take into account the secondary spectral dis- 
tortion in the radiation field, i.e. we cannot approximate the 
background radiation field J as a perfect blackbody spec- 
trum. This means, for example, that the extra photons from 
the recombination of Hel may reds hift into an energy range 
that can photoionize H( n = 1) IjDell' Antonio fc Rvbickil 
Il993t ISeaeer et al.ll200d) . We cab assess how significant 
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Figure 10. Comparison of the net 2p-ls (solid) and 2s-ls two- 
photon (dashed) transition rates of Hell as a function of redshift. 
The two-photon process is greater through most of the recom- 
bination epoch, so that most of the cosmological Helll — > Hell 
process happens through the two-photon transition. 




Figure 11. The normalized emission spectrum for the two- 
photon emission process (2 1 s— l^s) in Hel. Here y = v/v2 s —is, 
where i/ 2s -i s = 4.9849 X 10 15 H z. The crosses are the calculated 
points from iDrake et al ] ^196flft and lDrakel il98rt . while the line 
is a cubic spline fit. 

this effect might be by considering the ratio of the num- 
ber of CMB background photons with energy larger than 
-E 7 , n 7 (> to the number of baryons, tib, at different 

redshifts (see Fig. 1121 . 

Roughly speaking, the recombination of H occurs at the 
redshift when n 7 (> h p v a )/nB is about equal to 1. This is 
because at lower redshifts there are not enough high energy 
background photons to photo-ionize or excite electrons from 
the ground state to the upper states (even to n = 2), while 
at higher redshift, when such transitions are possible, there 
are huge numbers of photons able to ionize the n = 2 level. 
The solid line in Fig. 1121 shows the effect of the helium line 
distortions on the number of high energy photons (above 
Lycv) per baryon. The amount of extra distortion photons 
with redshifted energy larger than h p u a coming from the 
recombination of Hel is only about 1 per cent of the number 
of hydrogen atoms. Their effect is therefore expected to be 
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Figure 12. The ratio of number of CMB photons with energy 
larger than E 1 (n 7 (> Ej)) to number of baryons (?ib) is plot- 
ted against redshift z. The solid line includes the extra distortion 
photons from the recombination of Hel. From the graph, we can 
see that the recombination of H occurs approximately at the red- 
shift when the ratio of photons with energy > hu a to baryons is 
about unity. By the time the helium recombination photons are 
a significant distortion to the CMB tail above Ly a the density of 
the relevant photons has already fallen by 2 orders of magnitude, 
and so the effects can make only a small correction. 

negligibly small for x e . We neglect the effect of the helium 
recombination photons on the hydrogen line distortion, since 
it is clearly going to make a small correction (at much less 
than the 10 per cent level). 

As well as this particular approximation, there have 
been some other recent studies which have suggested that 
it may be necessary to make minor mo difications to t he re- 
comb ination calculations presented in ISeager et alJ <ll999l . 
2000). Although these proposed modifications would give 
only small changes to the recombination calculation, it is 
possible that they could have much more significant effects 
on the line amplitudes and shapes. Recent papers have de- 
scribed 3 separate potential effects. 

In the effective three-level model, iLeung. Chan fc Chul 

(2004) argued that the adiabatic index of the matter should 
change during the recombination process, as the ionized gas 
becomes ne utral, giving s light diffe rences in the recombina- 
tion history. fPubrovich fc Grachev! (|200rJ l have claimed that 
the two-photon rate between the lowest triplet state and the 
ground state and that between the upper singlet states and 
the ground s tate should not be ignored in the recombination 
of Hel. And IChluba fc Sunvaevl i2005f) suggested that one 
should also include stimulated emission from the 2s state of 
H, due to the low frequency photons in the CMB blackbody 
spectrum. Even if all of these effects are entirely completely 
correct, we find that the change to the amplitude of the main 
spectral distortion is much less than 10 per cent. We there- 
fore leave the detailed discussion of these and other possible 
modifications to a future paper. 

4.2 Possibility of detection 

There is no avoiding the fact that detecting these CMB spec- 
tral distortions will be difficult. There are 3 main challenges 
to overcome: (1) achieving the required raw sensitivity; (2) 
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removing the Galactic foreground emission; and (3) distin- 
guishing the signal from the CIB. 

Let us start with the first point. We can estimate the 
raw sensitivity achievable in existing or planned experiments 
(even although these instruments have not been designed for 
measuring the line distortion) . Since the relevant wavelength 
range is essentially impossible to observe from the ground, 
it will be necessary to go into space, or at least to a balloon- 
based mission. One existing experiment with sensitivity at 
relevant wavelengths is BLAST (Devlin et al. 2004), which 
has an array of bolometers operating at 250 /im on a balloon 
payload. The estimated sensitivity is 236 mJy in 1 s, for a 
30 arcsec FWHM beam, which corresponds to XI\ — 1.2 x 
10" 7 W m 2 sr 1 . Comparing with equation l|27|l for the peak 
intensity, it would take ~ 10 7 such detectors running for a 
year to detect the line distortion. The SPIRE instrument on 
Herschel will have a similar bolometer array, but with better 
beamsize. The estimated sensitivity of 2.5 mJy at 5a in 1 
hour for a 17.4 arcsec FWHM beamsize (Griffin, Swinyard 
& Vigroux 2001) corresponds to Xh = 4.4 x 10~ 8 W m~ 2 sr _1 
per detector for the la sensitivity in 1 second. So detection 
of the line would still require ~ 10 6 such detectors operating 
for a year. 

These experiments are limited by thermal emission from 
the instrument itself, and so a significant advance would 
come from cooling the telescope. This is one of the main 
design goals of the proposed SAFIR (Leisawitz 2004) and 
SPICA (Nakagawa et al. 2004) missions. One can imagine 
improvements of a factor ~ 100 for far-IR observations with 
a cooled mirror. This would put us in the regime where 
arrays of ~ 10 4 detectors (of a size currently being manu- 
factured for sub-mm instruments) could achieve the desired 
sensitivity. 

One could imagine an experiment designed to have 
enough spectroscopic resolution to track the shape of the 
expected line distortion. The minimum requirement here is 
rather modest, with only X/SX ~ 10 in at least 3 bands. 
An important issue will be calibration among the different 
wavelengths, so that the non-thermal shape can be confi- 
dently measured. To overcome this, one might consider the 
use of direct spectroscopic techniques rather than filtered or 
frequency-sensitive bolometers. 

Another way of quoting the required sensitivity is to say 
that any experiment which measures the recombination line 
distortion would have to measure the CIB spectrum with a 
precision of about 1 part in 10 5 , which is obviously a sig- 
nificant improvement over what can be currently achieved. 
A detection of the line distortion might therefore naturally 
come out of an extremely precise measurement of the CIB 
spectrum, which would also constrain other high frequency 
distortions to the CMB spectrum. 

Some of the desig n issues involved in such an experi- 
ment are discussed by iFixsen fc Mather! B2002T) . They de- 
scribe a future experiment for measuring deviations of the 
CMB spectrum from a perfect blackbody form, with an ac- 
curacy and precision of 1 part in 10 6 . This could provide 
upper limits on Bose-Einstein distortion fj, and Compton dis- 
tortion y parameters at the ~ 10~ 7 level (the current upper 
limits for y and n are 15 x 10 -6 and 9 x 10 -5 , respectively; 
IFixsen et all ll 996). The wavelength coverage they discuss 
is 2-120 cm -1 (about 80-5,000 /im), which extends to much 
longer wavelengths than necessary for measuring the line 



distortion. The beam-size would be large, similar to FIRAS, 
but the sensitivity achieved could easily be 100 times better. 
An experiment meant for detecting the line distortion would 
have to be another couple of orders of magnitude more sen- 
sitive still. 

Turning to the second of the major challenges, it will 
be necessary to detect this line in the presence of the strong 
emission from our Galaxy. At 100 /im the Galactic Plane can 
be as bright as ~ 10 3 MJysr _1 which is about a billion times 
brighter than the signal we are looking for! Of course the 
brightness falls dramatically as one moves away from the 
Plane, but the only way to confidently avoid the Galactic 
foreground is to measure it and remove it. So any experi- 
ment designed to detect the line distortion will need to cover 
some significant part of the sky, so that it will be possible 
to extrapolate to the cosmological background signal. The 
spectrum of the foreground emission is likely to be smoother 
than that of the line distortion, and it may be possible to 
use this fact to effectively remove it. However, it seems rea- 
sonable to imagine that the most efficient separation of the 
signals will involve a mixture of spatial and spectral infor- 
mation, as is done for CMB data (see e.g. Patanchon et 
al. 2004). 

In the language of spherical harmonics, the signal we 
are searching for is a monopole, with a dipole at the ~ 10 -3 
level and smaller angular scale fluctuations of even lower 
amplitude. Hence we would expect to be extrapolating the 
Galactic foreground signals so that we can measure the over- 
all DC level of the sky. This is made much more difficult by 
the presence of the CIB, which is also basically a monopole 
signal. Hence spatial information cannot be used to separate 
the line distortion from the CIB. The measurement of the 
line distortion is therefore made much more difficult by the 
unfortunate fact that the CIB is several orders of magnitude 
brighter - this is the third of the challenges in measuring the 
recombination lines. 

The shape of the CIB spectrum is currently not 
very well characterised. It was detected using data 
from the DIRBE and FIRAS experiments on the 
COBE satellite. Estimates for the background (XI x) are: 
9nWm _2 sr -1 at 60 ^m (Miville-Deschenes, Lagache & 
Puget 2002); 23nWm" 2 sr _1 at 100 ^m (Lagache, Haffiner 
& Reynolds 2000); 15nWm" 2 sr _1 at 140 ^m (Lagache et 
al. 1999; Hauser et al. 1998); and llnWm" 2 sr _1 at 240 fim 
(Lagache et al. 1999; Hauser et al. 1998). In each case the 
detections are only at the 3-5<r level, and the precise values 
vary between different prescriptions for data analysis (see 
also Schlegel, Finkbeiner & Davis 1998; Finkbeiner, Davis 
& Schlegel 2000; Hauser & Dwek 2001; Wright 2004). The 
short wavelength distortion of the CMB, interpretted as a 
measurement of the CIB (Puget et al. 1996) can be fit with 
a modified blackbody with temperature 18.5 K and emissiv- 
ity index 0.64 (although there is degeneracy between these 
parameters), which we plotted in Fig. |3] 

The CIB is thus believed to peak somewhere around 
100 /im, which is just about where we are expecting the re- 
combination line distortion. The accuracy with which the 
CIB spectrum is known will have to improve by about 5 or- 
ders of magnitude before the distortion will be detectable. 
Fortunately the spectral shape is expected to be significantly 
narrower than that of the CIB - the line widths are similar 
to the Sz/z ~ 0.1 for the last scattering surface thickness, as 
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opposed to 8\/\ ~ 1 for a modified blackbody shape (po- 
tentially even wider than this, given that the sources of the 
CIB come from a range of redshift Az ~ 1). 

One issue, however, is how smooth the CIB will be at 
the level of detail with which it will need to be probed. It 
may be that emission lines, absorption features, etc. could 
result in sufficiently narrow structure to obscure the recom- 
bination features. We are saved by 2 effects here: firstly the 
CIB averaged over a large solid angle patch is the sum of 
countless galaxies, and hence the individual spectral features 
will be smeared out; and secondly, the far-IR spectral en- 
ergy distributions of known galaxies do not seem to contain 
strong features of the sort which might mimic the recombi- 
nation distortion (see e.g. Lagache, Puget & Dole 2005). As 
we learn more about the detailed far-IR spectra of individual 
galaxies we will have a better idea of whether this places a 
fundamental limit on our ability to detect the recombination 
lines. 

Overall it would appear that the line distortion should 
be detectable in principle, but will be quite challenging in 
practice. 



5 CONCLUSIONS 

We have studied the spectral distortion to the CMB due 
to the Ly a and 2s-ls two-photon transition of H and the 
corresponding lines of Hel and Hell. Together these lines 
give a quite non-trivial shape to the overall distortion. The 
strength and shape of the line distortions are very sensitive 
to the details of the recombination processes in the atoms. 
Although the amplitude of the spectral line is much smaller 
than the Cosmic Infrared Background, the raw precision re- 
quired is within the grasp of current technology, and one can 
imagine designing an experiment to measure the non-trivial 
line shape which we have calculated. The basic detection of 
the existence of this spectral distortion would provide incon- 
trovertible proof that the Universe was once a hot plasma 
and its amplitude would give direct constraints on physics 
at the recombination epoch. 



6 ACKNOWLEDGMENTS 

This research was supported by the Natural Sciences and 
Engineering Research Council of Canada. S.S. is supported 
by the Carnegie Institution of Washington. S.S. thanks John 
Bahcall for stimulating discussions. We thank Argyro Tasit- 
siomi for helpying clarify Figure 1. 



REFERENCES 

Boschan P., Biltzinger P., 1998, A&A, 336, 1 
Burdyuzha V.V., Chekmezov A.N., 1994, AZh, 71, 341 
Burgin M.S., 2003, AZh , 80, 9, 771; English transl., 2003, 

ARep, 47, 709 
Chluba J., Sunyaev R. A., 2005, A&A, 446, 39 
Dell' Antonio LP, Rybicki G.B., 1993, in ASP Conf. Ser. 

51, Observational Cosmology, ed. G.Chincarini et al. (San 

Francisco: ASP), 548 



Devlin M.J. et al., 2004, Astronomical Structures and 

Mechanisms Technology, ed. A. Joseph & D. Lemke, Proc. 

SPIE, Vol. 5498, p. 42 
Derevianko A., Johnson W.R., 1997, Phys. Rev. A, 56, 1288 
Drake G.W.F., Victor G.A., Dalgamrno A., 1969, Phys. 

Rev., 180, 25 
Drake G.W.F., 1986, Phys. Rev. A, 34, 2871 
Dubrovich V.K., 1975, Pis'ma Astron. Zh., 1, 10, 3; English 

transl., 1975, Sov. Astr. Lett., 1, 196 
Dubrovich V.K., Stolyarov V.A., 1995, A&A, 302, 635 
Dubrovich V.K., Stolyarov V.A., 1997, Astrophys. Lett., 

23, 565 

Dubrovich V.K., Grachev S.I., 2005, Astronomy Letters, 
31, 359 

Fahr H.J., Loch R., 1991, A&A, 246, 1 
Finkbciner D.P., Davis M., Schlegel D.J., 2000, ApJ, 544, 
81 

Fixsen D.J., Cheng E.S., Gales J.M., Mather J.C., Shafer 
R.A., Wright E.L., 1996, ApJ, 473, 576 

Fixsen D.J., Dwek E., Mather J.C., Bennett C.L., Shafer 
R.A., 1998, ApJ, 508, 123 

Fixsen D.J., Mather J.C., 2002, ApJ, 581, 817 

Goldman S.P, 1989, Phys. Rev. A, 40, 1185 

Griffin M.J., Swinyard B.M., Vigroux L., 2001, The 
Promise of the Herschel Space Observatory, ed. GL. Pil- 
bratt, J. Cernicharo, A.M. Heras, T. Prusti & R. Harris, 
ESA Special Publication 460, 37 

Hauser M.G. et al., 1998, ApJ, 508, 25 

Hauser M.G., Dwek E., 2001, Annu. Rev. Astron. Astro- 
phys., 39, 249 

Hummer D.G., 1994, MNRAS, 268,109 

Kholupenko E.E., Ivanchik A.V., Varshalovich D.A., 2005, 
www.ioffe.ru/astro/QC/CMBR/csmn_nc.pdf 

Lagache G., Abergel A., Boulanger F., Desert F.X., Puget 
J.-L., 1999, A&A, 344, 322 

Lagache G., Haffner L.M., Reynolds R.J., Tufte S.L., 2000, 
A&A, 354, 247 

Lagache G., Puget J.-L., Dole H., 2005, ARA&A, 43, 31 

Leisawitz D., 2004, AdSpR, 34, 631 

Leung P.K., Chan C.W., Chu M.C., 2004, MNRAS, 349, 2, 
632 

Lipeles M., Novick R., Tolk N., 1965, Phys. Rev. Lett., 15, 
690 

Lyubarsky Y.E., Sunyaev R.A., 1983, A&A, 123, 171 
Martinis M., Stojic M., 2000, FIZIKA A (Zagreb), 93, 115 
Mather J.C., Fixsen D.J., Shafer R.A., Mosier C, Wilkin- 
son D.T., 1999, ApJ, 512, 511 
Miville-Deschenes M.-A., Lagache G., Puget J.-L., 2002 

A&A, 393, 749 
Nakagawa T. et al, 2004, AdSpR, 34, 3, 645 
Padmanabhan T., 1993, Structure Formation in the Uni- 
verse, Cambridge Univ. Press, Cambridge, UK 
Patanchon G., Cardoso J.-F., Delabrouille J., Vielva P., 

2005, MNRAS, 364, 1185 
Peebles P.J.E. 1968, ApJ, 153, 1 
Peebles P.J.E., Seager S., Hu W., 2000, ApJ, 539, LI 
Pequignot D., Petitjean P., Boisson C, 1991, A&A, 251, 
680 

Press W.H., Flannery B.P., Teukolsky S.A., Vetterling 
W.T., 1992, Numerical Recipes in C: The Art of Scien- 
tific Computing, Cambridge Univ. Press, Cambridge, UK 

Puget J.-L., Abergel A., Bernard J. -P., Boulanger F., Bur- 



12 W. Y. Wong, S. Seager and D. Scott 



ton W. B., Desert F. -X., Hartmann D., 1996, A&A, 308, 
L5 

Rubino-Martin J. A., Hernandez-Monteagudo C, Sunyaev 

R.A., 2005, A&A, 438, 461 
Santos J. P., Parente F., Indelicato P., 1998, Eur. Phys. J., 

D3, 43 

Schlegel D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 
525 

Scott D., Smoot G., 2004, The Review of Particle Physics, 

S. Eidelman et al., Physics Letters B. 592,1 
Seager S., Sasselov D.D., Scott D., 1999, ApJ, 523, LI 
Seager S., Sasselov D.D., Scott D., 2000, ApJS, 128, 407 
Spergel D.N. et al., 2003, ApJS, 148, 175 
Spitzer L., Greenstein J.L.Jr., 1951, ApJ, 114, 407 
Switzer E.R., Hirata CM., Phys. Rev. D, 72, 083002 
Tung J.H., Ye X.M., Salamo G.J., Chan F.T., 1984, Phys. 

Rev. A, 30, 1175 
Wright E.L., 2004, New Astronomy Reviews, 48, 465 
Zel'dovich Y.B., Kurt V.G., Sunyaev R.A., 1968, Zh. Eksp. 

Teor. Fiz., 55, 278; English transl., 1969, Soviet Phys. - 

JETP Lett., 28, 146 

This paper has been typeset from a TffjX/ DTjtX file prepared 
by the author. 



